library(here)
dir = here()
dir
wd = gsub("Figure", "Data/Case Studies/BGH", dir)
setw(wd)
setwd(wd)
d
wd
wd = gsub("Figures", "Data/Case Studies/BGH", dir)
setwd(wd)
d1 = fread("BGH_human_capital.csv")
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
setwd(dir)
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
View(d1)
exp(-1.37)
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
library(grid)
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
png(filename = "BGH_human_capita.png", height = 6, width = 7.5, unit = "in", res = 600)
grid.arrange(
rbind(
cbind(gA, gB, size = "first"),
cbind(gC, gD, size = "first"),
size = "first"
)
)
dev.off()
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
library(ggplot2)
library(gridExtra)
library(here)
library(data.table)
library(grid)
text.size = 10
dir = here()
wd = gsub("Figures", "Data/Case Studies/BGH", dir)
setwd(wd)
d1 = fread("BGH_human_capital.csv")
r = lm(d1$log.norm.pay ~ d1$edyrs)
r.square = summary(r)$r.squared
x = d1$log.norm.pay
y = d1$edyrs
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(10, 30 length.out = 50)
x = seq(10, 30, length.out = 50)
x = data.frame(x)
x = d1$log.norm.pay
y = d1$edyrs
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(10, 30, length.out = 50)
x = data.frame(x)
regression = predict(r, x,  interval = "prediction", level = 0.95)
View(regression)
x = d1$log.norm.pay
y = d1$edyrs
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(10, 30, length.out = 50)
x = data.frame(x)
regression = predict(r, x,  interval = "prediction", level = 0.95)
View(regression)
regression = data.frame(x = x, regression)
intensity.r2 =  paste("R^2 == ", round(r.square,2))
edu.plot = ggplot() +
geom_ribbon(data = regression, aes( x = x, ymin = lwr, ymax = upr), fill = "black", alpha = 0.05) +
geom_line(data = regression, aes(x = x, y = fit), col = "grey50") +
geom_jitter(data = d1, aes(x = edyrs, y= exp(log.norm.pay)), size = 0.3, col = "dodgerblue4") +
scale_x_continuous("Years of Education", breaks = seq(0, 40,2)) +
scale_y_log10("Normalized Pay", breaks = c(0.2, 0.5, 1, 2, 5, 10)) +
ggtitle("A.  Education") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5, vjust = -0.2),
axis.line = element_line(color = "black"),
axis.title.x=element_text(size=rel(0.9)),
axis.title.y=element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
edu.plot
x = d1$log.norm.pay
y = d1$edyrs
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(10, 30, length.out = 50)
x = data.frame(x)
regression = predict(r, x,  interval = "prediction", level = 0.95)
regression = data.frame(x = x, regression)
View(regression)
log(70)
log(17)
x = d1$log.norm.pay
y = d1$edyrs
r = lm(y ~ x)
predict(r)
plot(predict(r))
plot(x,y)
# regression
x = d1$edyrs
y = d1$log.norm.pay
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(10, 30, length.out = 50)
x = data.frame(x)
regression = predict(r, x,  interval = "prediction", level = 0.95)
regression = data.frame(x = x, regression)
intensity.r2 =  paste("R^2 == ", round(r.square,2))
edu.plot = ggplot() +
geom_ribbon(data = regression, aes( x = x, ymin = lwr, ymax = upr), fill = "black", alpha = 0.05) +
geom_line(data = regression, aes(x = x, y = fit), col = "grey50") +
geom_jitter(data = d1, aes(x = edyrs, y= exp(log.norm.pay)), size = 0.3, col = "dodgerblue4") +
scale_x_continuous("Years of Education", breaks = seq(0, 40,2)) +
scale_y_log10("Normalized Pay", breaks = c(0.2, 0.5, 1, 2, 5, 10)) +
ggtitle("A.  Education") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5, vjust = -0.2),
axis.line = element_line(color = "black"),
axis.title.x=element_text(size=rel(0.9)),
axis.title.y=element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
edu.plot
x = d1$edyrs
y = d1$log.norm.pay
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(10, 30, length.out = 50)
x = data.frame(x)
View(x)
d1 = fread("BGH_human_capital.csv")
# regression
x = d1$edyrs
y = d1$log.norm.pay
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(10, 30, length.out = 50)
x = data.frame(x)
regression = predict(r, x,  interval = "prediction", level = 0.95)
regression = data.frame(x = x, regression)
intensity.r2 =  paste("R^2 == ", round(r.square,2))
View(regression)
plot(x, y)
x = d1$edyrs
y = d1$log.norm.pay
plot(x, y)
# regression
x = d1$edyrs
y = d1$log.norm.pay
r = lm(y ~ x)
r.square = summary(r)$r.squared
x = seq(10, 30, length.out = 50)
x = data.frame(x)
regression = predict(r, x,  interval = "prediction", level = 0.95)
regression = data.frame(x = x, exp(regression))
intensity.r2 =  paste("R^2 == ", round(r.square,2))
edu.plot = ggplot() +
geom_ribbon(data = regression, aes( x = x, ymin = lwr, ymax = upr), fill = "black", alpha = 0.05) +
geom_line(data = regression, aes(x = x, y = fit), col = "grey50") +
geom_jitter(data = d1, aes(x = edyrs, y= exp(log.norm.pay)), size = 0.3, col = "dodgerblue4") +
scale_x_continuous("Years of Education", breaks = seq(0, 40,2)) +
scale_y_log10("Normalized Pay", breaks = c(0.2, 0.5, 1, 2, 5, 10)) +
ggtitle("A.  Education") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
legend.position = "none",
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5, vjust = -0.2),
axis.line = element_line(color = "black"),
axis.title.x=element_text(size=rel(0.9)),
axis.title.y=element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
edu.plot
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
library(ggplot2)
library(gridExtra)
library(here)
library(data.table)
library(grid)
text.size = 10
dir = here()
wd = gsub("Figures", "Data/Case Studies/BGH", dir)
setwd(wd)
d1 = fread("BGH_human_capital.csv")
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_rank.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_rank.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_rank.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_rank.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_rank.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_rank.R')
library(data.table)
library(magrittr)
library(ppcor)
d = fread("BGH.csv", stringsAsFactors = F)
d$Year = paste("19", d$year, sep = "")
d$Year = as.numeric(d$Year)
d$level = as.numeric(d$level)
# remove BGH outlier years
remove = c(1986, 1987, 1988)
d = subset(d, !Year %in% remove)
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_capital.R')
source('~/Desktop/human_capital/Supplementary Material/Figures/BGH_human_rank.R')
d1 = fread("Case Results.csv")
library(ggplot2)
library(gridExtra)
library(here)
library(data.table)
library(grid)
text.size = 10
dir = here()
wd = gsub("Figures", "Data/Case Studies", dir)
setwd(wd)
d1 = fread("Case Results.csv")
View(d1)
library(ineq)
d = fread("Case Results.csv")
d = fread("Case Results.csv")
View(d)
test = d[, ineq(Mean), by = .(Year, Source)]
View(test)
hist(test$V1)
View(d)
View(test)
plot(density(test$V1))
plot(density(test$V1, na.rm = T))
View(d)
d = d[Source == "BGH"]
test = d[, ineq(Mean), by = .(Year, Source)]
hist(test$V1)
plot(density(test$V1, na.rm = T))
View(d)
